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Abstract 

In some newly discovered materials the ratio of phonon to electron ener- 
gies is no longer small. We have investigated the basement of the recently 
proposed gauge-invariant self-consistent method and found conditions of its 
applicability. 

When describing electron-phonon interaction in usual metals, so-called adiabatic 
approximation is used. It assumes that due to electroneutrality fast electrons must 
follow "slow" ions. Narrow bandwidths and strong electron-phonon interaction make 
the orders of these velocities comparable. A question appears: Does this result in a 
breakdown of adiabaticity and increase of interelectron interaction (and in the result 
T c )? 

Traditionally such corrections were considered to be small as the Migdal param- 
eter Q ph /W was small |[ |3|, ||. But recently there were discovered materials 
(for example, fullerenes and high T c superconductors) where Q p h ~ W . So Migdal's 
theorem about smallness of adiabatic corrections |J] is violated. This gives the 
grounds for not just discussing these corrections, but even for considering an an- 
tiadiabatic limit Q p h 3> W. For example, it is stated that this results in increase 
of interelectron interaction and can increase T c for superconductors ||. But even 
with small Q p h/W there are contradictory conclusions about importance of such 
corrections (the author of Ref. |J considers them to be negligible, in contradiction 
to the authors of Refs. 0,0). 

There are two methods to take nonadiabatic corrections into account. One of 
them, which can be called Migdal's one, is based on the solution to the Bethe- 
Salpeter equation for the vertex function ^ |2|, |3|, [/J. In the lowest approximation, 
the first correction to the unit vertex is determined by the diagram in Fig. 1 (a). 
Then with the found vertex function the self-energy E is calculated (see Fig. 1 (b)). 
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This traditional method does not allow to take into account higher order diagrams 
analytically because of complexity of integration over momenta. 

Recently Y. Takada offered another method which he called the gauge- 
invariant self-consistent (GISC). In this method, based on the Ward identity 0, 
the vertex function is chosen as a functional of the self-energy which is supposed to 
be independent from momentum. So, when comparing this method with Migdal's 
one, we discuss validity of such approach. The parameter fl ph /W is considered to 
be small. 

For simplicity in the Letter we consider the model of the Fermi liquid with the 
usual Frohlich Hamiltonian M and the Einstein phonon spectrum. This corresponds 
to a 2 Fg;(fi) = ^XQphS(Q — Q p h) where Q p h is the phonon frequency independent 
from momentum, and A is the mass-enhancement parameter. It is assumed that the 
density of states is the constant iV(0) for — W < e < W where 2W = 2p F vp is the 
bandwidth. The chemical potential equals zero which corresponds to a half-filled 
band. 

We are interested in the quasiparticle self-energy E on the imaginary axis which 
is determined by the diagram in Fig. 1 (b). The phonon Green's function Do(iu u ) = 
— flp^cul + flpf^^ 1 , to v = 2vixT ', and the electron Green's function G(iu) n ,p) = {iuj n — 
Bp- — E(iuj n ,fi))~ l , uj n = (2n + l)vrT, is the electron spectrum. 

Assuming that E(zo; n , k) weakly depends on k, we take \k\ = pf = const, and so 
E = T,(iu n ). Then the diagram in Fig. 1 (b) corresponds to [[| 



E(iu n ) = I dVLa 2 F E {tt)TY, 



2Q 



- ul + W N(0) 

x G(iu n - iu) v , k - q)T(iu n , iu n - iu u , q) (1) 
q 

Takada's and Migdal's methods differ by a choice of the vertex function, and 
this is discussed below. We will obtain the results for T = and the variables will 
change as follows: uo n — > u, uo v — > uo' remaining on the imaginary axis. 

According to Takada's GISC method, the vertex function can be chosen as a 
functional of the self-energy Tt = IV [Ey]. This choice is based on the Ward 
identity [§ 

iuj u T(iuj n , iuj n — ioo u , k,k — q) — qT(iuj n , ioo n — iu) u , k,k — q) 

= G^iiuJn, k) - G^iiun - iuj u , k - q), (2) 

valid for all uj v and q. It is essential that the phonon Green's function does not 
depend on momentum. Only in this case one may use the Ward identity for finite 



momentum q 0. In Ref. || Takada proposes to use this identity for Vp\q\ -C 
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and choose || 

T T E = r T (^ n , - zuv) = 1 + ( — : — -), (3) 

lioo n 2[zu n — iu v ) 

which corresponds to the estimates of Refs. 0, |3| in this limit. 

The set of equations flU), @ is proposed to be solved by iterations. At the first 
step G = G (0) , T = = 1 give E T = E$! ] (Fig. 2 (a)). For T = 

E^P (za;) = —i\Q p h arctan — h iXQ p h arctan — — — — , (4) 

ilph W + Mph 



which corresponds to the results of Refs. 0, ]5|. At the second step G = G\Ejr], 
and the correction for Q p h/W <C 1 gives: 

Sy (zo;) = iX 2 Q p hb(m), where m = u/Q p h, b(m) = &i(m) + 62 (wi)- (5) 



bAm) = arctan 2 m, (6) 

2m 



7 / \ 1 /* 7 1 arctan! m — y) . . . . . 

62 (m) = - / dy sign(y-m), y = u U ph . (7) 

4 J y + 1 m — y 

—00 

So 



KmH-j^&hH^ . (8) 



4 m' 



Here Ej^ is obtained by using (H) with infinite limits of integration over de. 

(2) 

During this calculation the vertex is of the order of A for any q which is incorrect. 
So in Takada's method the leading term of the corrections does not have the small 
parameter Q p h/W, that is in this approach Migdal's theorem does not hold for the 
self-energy. This can be explained by the fact that the region vp\q\ <C \uj'\ where 
T ~ A gives in reality a small contribution to E. If, on the contrary, Vp\q\ ^> |u/|, 
the vertex T ~ XVt p h/W which gives E^ 2 ) ~ XVt^ h /W . At the same time, the sign 
of the vertex correction to the self-energy in GISC method is opposite to that in 
Migdal's one (see below). 

In the method which can be called Migdal's one, there is considered the first 
correction to the unit vertex function r^. It is shown on Fig. 1 (a). It has been 
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estimated in many papers [[I], [3] , but now it is becoming especially significant || as 
it has appeared that the parameter XQ p h/W can be not very small. 
The diagram in Fig. 1 (a) corresponds to the expression 



r(%- • • -a 7 \ 7,n ^ 2 F E (Q)2Q 

i ( . ^ - . Wv , q) = t E E j dn Nmuli + Q2) 



[iu n ,iu n - IU V 

xG^°\ioj n — iuj u i,p — q ')G^°\iuj n — iuj v i — iu u ,p — q' — g). (9) 



Assuming that | p — <f '| ~ and expanding e(p — q' — q) ~ — 9 ') — (J^f cos 
we get for T = 0: 

r^fza;, ic<j — ioj', q) = X— — arctan (arctan — arctan — — 

qv F uj' ilph il ph 

UJ uj — uj' 

-arctan— — h arctan — — — - ). (10) 

At Vp\q\ <C \uj'\ this expression corresponds to the Ward identity (Q) and to the 
choice of a functional (|3|). For Vp\q\ ^> \u)'\ the correction T^(iuj,iuj — iuj',q) ~ 
XQph/qvF ~ Afiph/W 7 , that is Migdal's theorem for the vertex function does hold in 
this method. 

The corresponding self-energy then equals (see Fig. 2) 

£W is equal to as it also corresponds to = 1. For the correction due to 
r( 2 ), S^ 2 ), one can get at Vt ph /W < 1: 

Q 2 uj 
H^(iu) = iX 2 '-rrrdJm) where m 



w Vv "' " "" n 



ph 



7T f 1 

dJm) = — / dy— sign (y) sign (y — m)(arctanm — arctan(m — y)). (11) 

4 J v + 1 



y* + l 
In the limiting cases 



4(m) = f 'r'.t • < i2 > 

1 m ^> 1 



4 m ' 



The sign of the correction £(, 2 ) is opposite to that of which favour tendency 
towards instability [0, f|. 
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In this approximation one can find the vector vertex function T and show that it 
really cannot be neglected. Using electron-hole symmetry, in our model, T(iu,iuj — 
iJ, k,k-q) = f« + f( 2 ) where f« 



m ' 



uj' Q^Fs, OJ UJ — uj\ 



qT^ = —i\Q P h(l arctan )(arctan — arctan 

qv F oj' il ph U ph 

The diagram in Fig. 1 (a) corresponds to if one substitutes by and 
I^ 2 ) (ioj n , icj n — ioj u , q) - by (ico, iuj — iuj' , —q) . 

If one designates x = u'/qvp, then using at Q ph /W <C 1 the expression flTU| ) 



for T( 2 \ we can find that the terms wT" and qT^ from (^|) become of the same 
magnitude if xq arctan(l/xo) = 1 — xq arctan(l/xo), from which x = 0.43. If, 
however, u'/qvF < Xq, the vector term dominates and it may not be neglected. 

So in the Letter there is consistently considered contribution of nonadiabatic 
effects to the electron self-energy in first orders in A. The two methods of taking 
into account the nonadiabatic corrections were compared: 1) Migdal's and 2) the 
GISC proposed by Takada. 

It is shown that the latter gives for the self-energy physically incorrect result. 
This is connected with the fact that the leading contribution to S is made by the 
region qvp 3> to' where V ~ XQ p h/W. This means that in the Ward identity one may 
not neglect the vector term as proposed by Takada. In the frame of the standard 
(Migdal's) approach, in the lowest order, there is evaluated the parameter x = 
uj'/qvF at which contributions of the vector and scalar terms become the same (see 
Fig. 3). They are equal at xo = 0.43. As it was shown above, the region of large 
q's indeed gives the main contribution to the self-energy. That is why the GISC 
method, which works only for small q's, may not be applied. 
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by ISI Foundation and EU INTAS Network 1010-CT 930055. The authors are also 
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Figure Captures 

Fig. 1. (a) The equation for the vertex function r^ 2 ^ in Migdal's method. The 
outgoing lines are shown for clearness and are not included in the definition of T. (b) 
The equation for the self-energy S. The wavy line is a bare phonon Green's function 
D , the double solid line is an electron Green's function G, and the hatched triangle 
is a vertex function T. 

Fig. 2 (a), (b). Diagrams which are included when calculating the self-energy 
S. The solid straight line is a bare electron Green's function G^. 

Fig. 3. The ratio of the vector to the scalar terms in the Ward identity. The 
GISC method may only be used for large x's (small q's). 




S(iw n ) = 




r (iu; n , iu n - ita 



iw n -iu„, k- q 



(b) 



\u w q 




